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We have considered gravity in a five-dimensional warped product space-time, with a time- 
dependent warp factor and a time-dependent extra dimension. The field equations have been solved 
for a spatially flat FRW brane and the energy conditions and the nature of bulk geometry have been 
examined. At low energies, the trapping of fields within the brane implies a correlation between the 
warp factor and the extra-dimensional scale factor. Generally, the bulk is not conformally flat. At 
high energies, the bulk is sourced by a dilaton-like minimally coupled scalar field, which depends on 
the warp factor. Both the bulk and the brane are found to be Anti de Sitter, with some exceptions. 



I. INTRODUCTION 



The theories of particle physics beyond the standard model |l|-|4j , taught us that the universe is possibly a (3+1)- 
dimensional spacetime embedded in a higher-dimensional bulk, with the extra-dimensions playing a crucial role in 
a dynamical spacetime. The knowledge that the gauge fields may consistently appear on the ten-dimensional (10D) 
boundary of a .^-symmetric eleven-dimensional (11D) space-time @ and that the effective five-dimensional (5D) 
theory of the strongly coupled heterotic strings is a gauged version of the (N — 1) 5D supergravity with four- 
dimensional (4D) boundaries jg], have led to extensive studies in this area. The simplest illustration can be given in 
the framework of a 5D theory with a single warped extra dimension 0-Q ■ The issue of localization has been addressed 
by many [Iol[I3 |. Cosmological consequences have been studied in the M-theory models with two branes at finite 
distances apart, acting as the boundaries of the bulk spacetime [HHHl, anc ^ m tne Randall- Sundrum (RS) model [l(J 
of both the single brane and two brane versions [15l - [l9| . 

The bulk is "slightly warped" in the scenarios of heterotic M theory with the inter-brane distances being represented 
. by scalar fields (called moduli), introduced through the integration constants in the classical solutions. For solutions 
existing over a continuous range of integration constants, the corresponding moduli are massless, giving rise to long 
range scalar forces, which tend to evolve cosmologically, being time-dependent. This leads to time-variations of the 
fundamental constants [2fj| . indicating that at least some moduli fields are not stabilized. To avoid this, the moduli 
fields are assumed to be stabilized by some mechanism in the early universe 21]. On the other hand, the RS model is 
"highly warped" and without any bulk scalar (except for the cosmological constant), with matter fields on the brane 
being responsible for the dynamics of the brane. However, string theory requires the dimensionally reduced low energy 
effective action in 5D, to include both the higher-order curvature terms, as well as dilatonic scalar fields Jc, [22 1. Thus 



scalar fields were introduced in the bulk of the RS model to stabilize the scale of the extra dimension [23j , with the 
single modulus, called the "radion", being related to the thickness of the AdS slice. The resulting 4D theories at low 
energies in both the M theory and RS models are very similar to the scalar- tensor theories [IH, [2J, [25[ . 

In the RSI setup without any bulk scalar, the weak gravity at low energy on either brane is Brans-Dicke type. The 
inter-brane distance is arbitrary and the corresponding mode is massless, thereby playing the role of a Brans-Dicke 
scalar. In presence of a bulk scalar field which couples to the branes, the inter-brane distance is no longer arbitrary 
and the expected massless mode acquires a nonzero mass, so that Einstein gravity is recovered [13, [26l. In the RS2 
setup, although, the hierarchy problem cannot be solved, there are several important consequences (27H29J. One of 
them is an alternative scenario of the inflationary universe, where slow-roll inflation on the brane is driven solely by 
the dynamics of a dilaton-like minimally coupled scalar field living in the 5D bulk, and not by an inflaton field. The 
5D effective action includes scalar fields of gravitational origin, having a suitable potential and Newtonian gravity is 
recovered in the 4D effective theory. The bulk scalar also facilitates to generate the hypersurface in the form of the 
so-called thick brane [3fj| . 

In the RS model, the bulk and the branes are maximally symmetric. The fields living on the negative tension 
brane, have masses proportional to the radion expectation value and for generic matter in the bulk, the one loop 
effective potential is not sufficient to stabilize the radion at a phenomenologically acceptable mass [3lj . In the warped 
braneworld models which are not maximally symmetric, the bulk contains a dilaton field, acting as the inflaton field, 
with an exponential potential coupled to the 5D gravity. The solutions have a power-law warp factor a(y) ex y q , where 
y is the distance along the extra dimension and q is determined by the parameters in the Lagrangian. This yields a 
power-law inflation on the brane [32|, |33j , although the bulk itself is not inflating. The RS model is recovered in the 
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limit q — ¥ oo. The back reaction of the dilaton coupling in the bulk to the 4D brane gives rise to the non- vanishing 



component of the energy- momentum tensor along the fifth dimension 34| . 

Many authors have constructed various generalized and time-dependent solutions in an attempt to develop realistic 
models which can address the unsolved issues in astrophysics and cosmology (35l-l4l| . Cosmological compactifications 
on Einstein spaces of negative curvature, have yielded accelerating 4D FLRW cosmologies [331. The no-go theorem 
forbidding acceleration in the time- independent compactifications of string theory, are not applicable to a general class 
of warped supergravity models that allows at least one noncompact extra dimension [42[ . The de Sitter solutions of 
warped supergravity have finite 4D Newton's constant and the 4D massless graviton wave function is normalizable on 
an inflating brane (FLRW universe). Recently, dynamical brane solutions in supergravity and in string theory have 
also been studied [43|]. Depending on the ansatz for the fields, the warp factors in the solutions can be time-dependent. 
After compactification of the internal space, such solutions lead to the FLRW cosmology with power-law expansion. 
However, the power of the scale factor is too small for the solutions to yield a realistic expansion law. This opens up 
the possibility of further investigations with generalized spacetime geometries. 

With this aim, we have considered a thick brane [3] embedded in an AdS bulk, in presence of a dilaton-like bulk 
gravitational scalar which is time-varying, a general exponential warp factor which depends both on time as well as 
on the extra-dimensional coordinate (a special type was considered in [39]) and a time-dependent extra dimension. 
During the early stages of cosmological evolution (the high energy regime) , the universe underwent a very rapid phase 
of expansion, so that the process of localization of gravity was presumably time-dependent. This can be described by 
a time-dependent warp factor, which in turn gives rise to a time-varying bulk cosmological constant, with the scale of 
the extra dimension also time-dependent. The gravitational force law was sensitive to the background cosmological 
expansion [27| . and the 4D Newton's constant was time- varying. Gravity, along with scalar fields and particles 



could access the extra dimensions [45[. The effective cosmological constant of the hypersurface is determined by 



the curvature of the bulk metric, which is time-varying due to the time-dependent process of localization, and this 
curvature controls the 4D Newton's constant. The dynamical evolution of the brane is determined by the Einstein 
equations for the combined brane-bulk system and the dynamics of the gravitational field is modified by the presence 
of the time-dependent warp factor. If this modification takes place at today's Hubble scale, Hq, then it is expected 
that this will alter the gravitational force law at distance scales much smaller than H^ 1 [46| . leading to an expansion 
history that can be identically reproduced by a dynamical dark energy model. Somewhere around the end of inflation, 
the time-dependent process of localization came to an end, leading to the finite 4D Newton's constant in the present 
universe. This can be related to the infinite extent of the fifth dimension and is possible if the universe is inflating. 
The warp factor is regular at the boundary of the 3-brane and is well-behaved at infinite distance from the brane. The 
brane tension is induced both by the bulk cosmological constant, as well as the curvature related to the expansion of 
the (3+1) spacetime and the physical universe eventually accelerates. For an AdS bulk, the zero- mode graviton may 
not be localized on the brane [47| . 

Having presented the mathematical preliminaries in Section II, we construct the five-dimensional field equations in 
Section III and examine the status of the different energy conditions in the bulk. To get a clear picture of the geometry 
and the corresponding physical implications, we have considered a number of toy models. We have considered an 
isotropic bulk and have determined the nature of the bulk matter and that of the warp factor. To determine the exact 
analytic solutions, we have assumed a definite correlation between the warping function and the extra-dimensional 
scale factor, following the results of supersymmetric brane models 48]. In Section IV, we assumed a vanishing flux 
term for the generalized metric, signifying confinement of matter or energy (or both) within the 4D hypersurface, 
and have investigated the possibility of having the bulk as a conformally flat spacetime and the requirements for a 
stabilized bulk. In Section V, we consider the flux term to be non-zero, which therefore represents the condition 
when gravity becomes essentially higher dimensional. In this regime, we assume that exotic matter in the form of a 
dilaton scalar (minimally coupled to gravity) propagates in the bulk along with gravitons. We have considered two 
types of scalar fields: a massless non-self-interacting scalar field, for which the solutions are obtained for various types 
of warping function, and a scalar field with self-interaction potential for the most general warped metric in 5D, for 
which the solutions are obtained and the nature of cosmology in the bulk and the brane are analysed. In Section 
VI, we have examined the effect of placing the brane within the bulk and have discussed the nature of localization of 
gravitons and other matter fields, in brief. The summary and discussions have been presented in Section VII. 



II. MATHEMATICAL PRELIMINARIES 



Let us consider the generalized 5D action in presence of a bulk scalar field with potential, minimally coupled to 
gravity, as follows [H El S H3 



S = / d 5 x^ 



_L(j? _ 2A) + \g AB V A ^ B ip - Vty) 



(1) 
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supplemented with the brane curvature term 




where g~AB is the 5D metric of signature (H -), A is the bulk cosmological constant and R is the 5D scalar 

curvature. The first term in (JTJ) corresponds to the 5D Einstein-Hilbert action in presence of the bulk scalar. The 
Lagrangian density L m represents all other contribution to the action which are not strictly gravitational, including 
the contribution of the matter fields localized on the brane and any interaction between the brane and the bulk. The 
constant R is related to the 5D Newton's constant Gr$) and the 5D reduced Planck mass M( 5 ) by the relation 

R 2 = 8^G (5) = M^. (3) 

The dynamics of the 5D space-time, including that of the 4D hypersurface representing the observed universe, is 
determined by the 5D field equations 

Gab = -AgAB + k 2 T A b (4) 

where Gab is the 5D Einstein tensor and Tab represents the 5D energy-momentum tensor. For a generalized 5D 
metric ansatz of the type as in [15| , there may exist definite correlations between the various scale factors in the line 
element. A specific correlation (viz. n 2 = a 2 = e 2 ?) leads us to the so-called "warped spacetime" [5l| . Study of 
such spacetimes is important in theories of gravity, since they constitute a wide variety of exact solutions to the field 
equations. Such warped spacetimes has been used to solve the hierarchy problem and to recover Newtonian gravity 
in the low energy effective theory • We therefore choose a 5D metric ansatz given by 

dS 2 = e 2f(t,y) (^2 _ R 2 {t){dr 2 + r 2 dd 2 + ^^2^2)) _ ^ y)dy 2 (g) 

where y is the coordinate of the fifth dimension, t denotes the conformal time and the function A(t, y) parametrizes 
the scale of the extra dimension at different times and at different locations in the bulk. The "warping function", 
/, is a smooth function and e 2 ^ t,y ^ is the time-dependent warp factor. The observed universe is represented by the 
hypersurface at y = ya. The complicated nature of the field equations owing to the generalized metric makes it 
necessary to impose additional constraints on the geometry to determine the solutions. 

The RS scenario led to unnatural correlation between the bulk cosmological constant and the brane tensions. 
However, brane-bulk supersymmetry correctly correlates the brane tensions and the bulk cosmological constant in 
the supersymmetric RS models [48|. Under this condition, the metric §5$ will be a solution of the 5D Einstein 
equations, with the warp factor being related to the extra- dimensional scale factor. For a stabilized bulk, A{t, y) must 
settle down to a fixed value. A stabilized bulk with constant curvature and characterized by a negative cosmological 
constant is free from the appearance of nonconventional cosmologies [52|. For a conformally flat bulk (a space of 
constant curvature), the extra-dimensional scale factor can be a growing function of time [22]]. For such a bulk, the 
resulting field equations are much simpler owing to the simplified geometry. However, a generalized bulk may not be 
conformally flat. A spatially flat FRW-type braneworld is appropriate for achieving cosmological solutions and can be 
embedded in any constant curvature bulk (53[, in spite of representing a dynamical brane. In our model, we extend 
it to the case of any arbitrary bulk which contains the present universe embedded in it and evolving with time. 



III. FIELD EQUATIONS AND ENERGY CONDITIONS IN THE BULK 

The non- vanishing components of the 5D Einstein tensor for the space-time under consideration are 
& _ 3 (R 2 2Rf fA RA\ 3 ( 2fl2+f , f'A'\ 



(7) 



(8) 
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6^2 

A 2 : 



(9) 



and 




ARf fA 




f'A' 



(10) 



Above, a dot represents differentiation with respect to the conformal time t and a prime stands for differentiation 
with respect to the fifth coordinate y and I = r,6, 4>. The components of the bulk stress-energy tensor are given by 



ft = p, 



rpi _ _p 



TV — —P 

v ~ v> 



T$ = Q. 



where p, P and P y are the energy density and pressure in the bulk, with Q representing a non-trivial flux term. 
The non-zero G* term with the corresponding T* component, indicates a flux of matter or energy (or both) from 
the 4D hypersurface into the bulk along the extra dimension. The energy scale necessary to ensure such an access 
is few hundred GeV. For such a generic case of the bulk, the corresponding energy-momentum tensor may not be 
diagonalizable, although for "physically reasonable" matter content on the hypersurface [54| (and may be in the bulk), 
we expect the stress-energy tensor to be diagonalizable. In the previous work (39[, it was found that the validity of 
the strong energy condition (SEC) in the bulk is governed by the nature of warping, as well as by the effect of the 
extra dimension. Keeping in mind the possibility of a negative pressure component of matter-energy (dark energy 
or cosmological constant), we require that at least the weak energy condition (WEC) and the null energy condition 
(NEC), must be obeyed. This implies that we must have 



p>0 and p + P t >0 (i = 1,2,3,4). 

Assuming 8nG^ — 1 in the remaining analysis, we obtain the field equations in 5D as 

-A + p, 



(11) 



(12) 



&j = -A-P, 



(13) 



and 



G„ 



-A - P„ 



-A - Q. 



(14) 



(15) 



The complicated nature of the various components of the Einstein tensor makes it difficult for us to comment directly 
on the status of the energy conditions from the above equations. To proceed further, let us write down the energy 
conditions in terms of the energy eigen values (say, a^), which must be real for matter-energy obeying the WEC [551 ] . 
These eigen values are given by the roots of the characteristic equation for the energy-momentum tensor: 



This implies that 



det(T - at) = 0. 



-p + a Q 

P + a _ 

P + a _ 

P+a 



(16) 











P, 



0. 



Solving this equation, we find that the eigen values are given by 



a 0A 



(p-p y )±j(p+p v r-4Q 2 



and 



a\ = d2 = d3 = P. 



(17) 
(18) 
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For these eigen values to be real, we must have, 

(p + P y ) 2 -AQ 2 >0, (19) 
along with the following restriction, so as to satisfy (jlip . i.e. 

a >0 and a - a, t > (i = 1,2,3,4). (20) 

The conservation of the energy-momentum tensor T b a ; a = leads us to the two equations 

P + KP + P) (y + l) +(P + Pv)j=0, (21) 

and 

P y + (p-3P + 4P y )f = 0. (22) 
where = H, is the Hubble parameter on the brane. 

A. Solutions for an isotropic bulk 

For an isotropic bulk, P = P y , so that the limiting case for condition (fT9)) reduces to 

(p + P) = \2Q\, (23) 

where Q may be positive or negative but (p + P) must be non- negative so as not to violate the WEC. Assuming 
positive flux for Q, we can rewrite equation (|2T1) in the form 

P=-2g(^3/+^ + ^, (24) 

which represents the equation of continuity for this isotropic bulk. For a surface enclosing a given volume of the bulk 
(such that it encloses the 4D universe in the form of a 3-brane), since Q is positive, the energy density inside this 
volume should decrease with time, if the quantity inside the brackets is positive. Thus, the rate of decrease of energy 
density also depends on the rate of variation of the warp factor, the rate of expansion of the observed universe and 
the normalized rate of variation of the extra-dimensional scale factor. From (|22[) we get 

(P + P) = ~J- (25) 
The l.h.s. of the above equation will be non-negative under the following possibilities: 

i P' = 0, which implies that the pressure is a constant. This is possible when the bulk matter energy-momentum 

tensor is that of a perfect fluid representing dust. We can have /' to be either positive or negative, which means 
that the warp factor may be either growing or of decaying type. However, the decaying type is more favorable 
for cosmologically relevant solutions 39] . 

ii P' > 0, for which /' should necessarily be negative i.e. the warp factor is of decaying type, so that the WEC is 

valid. 

iii P 1 < 0, in which case /' is positive and the warp factor is of growing type. 
Combining equations (fT3"|) and (IT41 for this isotropic bulk we obtain the relation 

1 (R 2R 2 5Rf - 2 •• fA 2RA A\ 3 / d ( f 



V Z + f- V-^TT-T =-T IT. Pt • (26) 



e 2 f \R R 2 R J J A R A A A \dy \A 



Since the l.h.s. involves only time-derivatives and the r.h.s. only the derivative with respect to the fifth coordinate, 
without any loss of generality, we can assume that they are separately equal to some constant (say, C). From the 
r.h.s. we can infer that the variation of the extra-dimensional scale factor along the direction of the extra coordinate 
at different locations in the bulk, depends on the variation of the warp factor along the extra dimension. However, 
accurate predictions can only be made if we know the exact functional relations. 

Equating the constant C to zero, we obtain the following relation from the l.h.s. of (|2l)|) 



/ 2/--r \+f-—--- =0. (27) 
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B. A specific relation between / and A in case of the isotropic bulk 

Remembering that the warp factor may be related to the extra-dimensional scale factor, we assume that 

2f=j=B(t,y). (28) 
In that case the relation between / and A may be of the type / = ^InA + constant. 

1. A de Sitter brane 
For a dS 3-brane, the scale factor R(t) is given by 

R(t) ~ e m . (29) 
The l.h.s. of (|26|) can now be simplified further to yield: 

H + 'SH 2 + Hf + f- j=0, (30) 

and finally we have 

2H + 6H 2 + H B - 2B 2 - B = 0. (31) 

A straightforward solution of (J33J) is obtained when B is simply a function of time. In such a case, the spacetime is 
unwarped, with the extra-dimensional scale factor depending only on time. However, such type of solutions destroy 
the generalized feature of the spacetime geometry that we intend to investigate. The other way round is to try with 
various explicit functional form of B and see which one yields a physically consistent result. 

2. A radiative brane 

The scale factor for a radiative brane is 



Using ([28]t. we have from ([22 



R(t) ~ Vt. (32) 



Substituting (1521 in (1551 . we get 



/ + 4/ 2 1 



(34) 



1 + 2tf 4i 2 ' 

The above gives the nature of the time- variation of the warping function at a given location in the bulk, specified 
by y = constant. For a static brane, y = constant is sufficient to specify the position of the brane along the extra 
dimension. However, for a dynamic brane, the deviation of the hypersurface from the tangent plane will change with 
time [39l [53(1 . This in turn will monitor the nature of the localization of fields on the brane. Since the warp factor is 
smooth, we assume that the deviation of the hypersurface changes slowly with time. During the epoch when, / = 0, 
(|34p may be simplified further to yield 

/=f, (35) 

where a = (2 ± \/68)/32, leading to a warping function of the type 

/ - ln(t a cosh{cy)), (36) 

with the extra-dimensional scale factor as 

A - t a cosh(cy). (37) 

This is just a particular solution, for which the warp factor is of product type. The general solution when / ^ and 
with the deviation varying rapidly, will be more complicated. 
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IV. LOW ENERGY GRAVITY 



We know that at low energies, gravity is localized at the brane along with particles. To prevent matter or energy 
flowing out of the brane along the fifth dimension, we require the flux term to be zero, i.e. T* = 0, which implies that 
Gy = and hence we obtain 

./' = -/'■ (38) 

Assuming that both A, f and their first order derivatives are continuous, (|1 T[) can be easily integrated to give the 
result 

A(t,y)= X (y)f'(t,y). (39) 

Thus the extra- dimensional scale factor is indeed related to the warp factor. The extra-dimensional scale factor at 
a given y and t, depends on the way the warping function varies along the extra dimension at that instant at the 
given location and will be different at different locations in the bulk. The manner in which gravity is localized on 
the 4-dimensional hypersurfaces at different locations in the bulk (being monitored by the warp factor), will also be 
different at different times. 

Consequently, the Einstein equations for the above spacetime get reduced to the form 





and 



(41) 



p *+^ = -w(R + w + nr +f +f + 7- (42) 



R , R 2 3Rf i2 , z\ 6 

2 

x 



Thus we have six unknowns (/, i?, x, Pi P and P y ) and three independent equations. To proceed further, we need 
to impose additional constraints to solve the field equations. 



A. Possibility of a conformally flat bulk 

The nonzero components of the Weyl tensor for the spacetime under the condition A(t, y) — x(y)f '{t-> V)t are given 

by 

Cjg - «— x ( ^ (| - f (/ + |) + / - ? + £ (| + V) - f ) ) (43) 

where, only the constant factor varies for the different components. For this bulk to be conformally flat, we need the 
following condition to be satisfied: 

f(t,y) = F 1 (t)F 2 (y) (44) 
where F\ and F% are arbitrary functions such that 

dF 2 (y) , N 

y = 0. 45 
ay 

This implies that the conformally flat solutions have f(t, y) — (Fi(t) x constant) and hence A(t, y) — x{y)f'{^i v) = 0- 
Therefore, we cannot have a conformally flat bulk with such a generalised scale factor for the extra dimension, although 
we obtain the conformall y fl at 5D bulk if the extra-dimensional scale factor is given by A = A(t), which may be a 
growing function of time [33j , corresponding to a moving brane within a static bulk. 
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1. A particular case 



We now examine a very special type of warped metric for which f(t, y) = ln(at + bcosh(cy)), such that 

e 2f(t,y) = A 2 (t, y) = (at + bcosh(cy)) 2 . (46) 

In this case, the Einstein tensors for the bulk and its 4D part are diagonal, and the bulk turns out to be conformally 
flat with identical components for the Riemann tensor if the following condition is valid: 

|-|=0. (47, 

Under this condition, the Hubble parameter of the brane turns out to be a constant. The Weyl tensor of the 4D part 
vanishes identically. The Ricciscalar for the 5D bulk, and for the corresponding 4D part are identical and is given 
below: 

6 (& R 2 \ 

~~{at + bcosh(cy)f \7 + B? J ' ( ' 

indicating that both the 4D part and the bulk are AdS, provided > The bulk pressure is isotropic along the 
three spatial directions on the brane, but is different along the extra dimension. The 4D part itself is isotropic. The 
scale factor for the 3-brane is then given by R(t) = ae Ht , and its deceleration parameter is given by the unique value 

q = -l (49) 

representing a universe in an epoch of uniform acceleration. 



B. The case of a stabilized bulk 



For this condition, we must have A — 0, so that A — A(y). Re-scaling A(y) — > 1, we find that for T* = 0, we need 
/' = 0. Thus either / = f(t) or / - f(y). 



1- / = /(»)•■ 

This corresponds to the usual RS scenario. For /" = 0, we have / = const x y, which represents a dS 3-brane 
embedded within a 5D AdS. When the bulk is of constant curvature, the bulk Weyl tensor vanishes and the bulk is 
characterized by a negative cosmological constant with ps = —Pb- 



2- / = /(*): 



This is the unwarped case. We can absorb f(t) inside the 4D metric. 



3. f(t,y) = l(lnr(t) + lnT(y)) : 

This is a particular case where the warp factor is taken in the product form and has been discussed earlier (39| . The 
assumption was made in order to transform the Einstein equations into a simpler form, which can be solved by some 
straightforward calculations. The energy-momentum tensor has only the component other than the components 
on the brane. As mentioned before, the validity of SEC in the bulk for this metric is governed by the nature of 
warping and the effect of the extra dimension. For an isotropic bulk, with r(t) = t -1 / 2 and R(t) — t 1 / 2 , the WEC is 
valid. The 4D universe is initially decelerated, but makes a transition to an accelerated phase at later times, leading 
to an interpretation of dynamical dark energy. 
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V. GRAVITY AT HIGH ENERGIES 

Here the energy scale is in the GeV range and particles can escape into the extra dimension. We shall consider two 
specific cases of the bulk scalar field, namely an ordinary, non-self-interacting massless scalar field and a scalar field 
with potential. 

A. A massless scalar field in the bulk 

The energy-momentum tensor for a non-self-interacting massless scalar field source ip(t,y) [4l[ is given by 

Tff*« r = -d^dji, + \g u d K ^d K 4, (50) 

Case 1: 

Here we assume the warping function f(t,y) to be given by 

f(t,y) = ^(]nr(t)+lnT(y)). (51) 

With this choice, the 5D metric becomes 

dS 2 = T(t)T(y) (dt 2 - R 2 (t){dr 2 + r 2 d9 2 + r 2 sin 2 9d<p 2 )) - A 2 (t, y)dy 2 . (52) 

The difference from [39] is in the functional form of the extra-dimensional scale factor. The components of the Einstein 
tensor for the bulk are now 

- _ 3g tt ( ARr AR 2 r 2 2fi\ 6g tt (V" T'A'\ 

G ty =^ (54) 



- oy m[ 2r t 2 6Rf AR 2 AR\ 6g yy (T'\ 2 



r, _ 9ij 1 4r 3f 2 8Rf AR 2 8R 2tA 8RA AA\ 6g tJ (T" _ V A' . 
° y ~i^\7" ~ + _ R7 + "R^ + _ R + TT + lM + T y l AA^\Y ~TA 1 

The components of the bulk energy-momentum tensor for the given metric are 

f ' = -(^ + S) = -^ < 57 ' 

1 M *" 



and 



T l = ib^ - — — 

1 2tT V 2 A 2 



v tT ■ 



To obtain the solutions for the scalar field, let us assume that 

fl>(t,y)=M*) + Mv)- (58) 
From the iy-component of the field equations, we obtain 
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i'i=k^ and ip' 2 = -j^jr, 

where k is an arbitrary constant. Consequently, the solution for the bulk scalar is given by 



, . . , f 1 OA , 3 , 
ip(t,y) = k / ——— at — —mr + constant. 



(59) 



As a toy model if we assume that 

A(t,y) = A 1 (t)xA 2 (y), (60) 

then the result is 

3 

ip(t,y) = klnAi -InT + constant. (61) 

2k 



However, we must remember that the assumption (|60[) is made only to obtain a clear picture of the nature of solution 
for a specific case, at the expense of the general character of the spacetime. 

Equating the coefficient of -4j on both sides of the j/y-component of the field equations, we find that for a non-zero 

value of y~ we require 

r(y) = C ie ±c y. (62) 

Considering a radiative brane with r as a power-law function, for which 

R(t) ~ Vi, (63) 

and T{t)~t m , (64) 

we obtain 

Ai(t) - H ro+1)/2 , (65) 

A 2 (y)^l. (66) 

To obtain a physically meaningful solution, A\(t) must be a decreasing function of time, which is true if we choose 
m = — i. This choice of m is based on the results obtained earlier (39[. 

Case 2: 

Let us now consider f(t, y) = 2bty, where b is a dimensional parameter. In this case we have 

Gty =^ b (^J- ^ (67) 

and 

f t y = -fa'. (68) 

Once again assuming a solution of the type given by (|58j). we obtain the following relations for the scalar field ip: 

i'i = -h - l) and ^2 = f , 

for which, the solution is of the type 

, , , , / f t OA , \ 3by 

tp(t,y) — —ki\ / —— — at — t ] + —, h constant, (69) 

\J A dt J ki 

the exact solution being obtainable only when A(t, y) is explicitly defined. 
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B. A scalar field with potential in the bulk 

Let us now assume that the 5D bulk includes a dilaton-like background scalar field ip(t,y) with self-interaction 
scalar potential V{ip) (22l. [37|. The energy-momentum tensor for a minimally coupled scalar field with potential is of 
the type 



f scalar = _Q^ dj ^ + gjJ ( -g KL d K ^d L ^ + V ty) ) . (70) 



The equation for the bulk scalar field is given by 



D (5)V> + ^=0. (71) 
With the assumed bulk metric (0, the scalar field equation of motion reduces to the following form: 

- >" + c ' 2 ' ( 2/ + f + i) * + i (-»' + £) *' = ~w < 72 > 

while the Einstein equations in the bulk are 

e 2 f\R 2 + R 1 A RA A 2 { 1 +/ 



_3_ / R R?_ 3R£ ^ - 2 
\ R + R 2 +_ R~ 4 ' 



1 j 2R i R 2 ARf i - 2 ■■ fA 2RA A\ 3 , .„ 



7) - 


/ y- 2 

\2e 2 f 


V;' 2 
+ 2I 2 


-4'- 






j A 2 


y> 2 

~ 2e 2 / + 


2A 2 + 


- 7) 


y> 2 

~ 2e 2 / 


y/ 2 

2A 2 4 



// -5 //= T' (74) 

+ /" + /VS£ = £ + £ + V + A, (75) 



Equation (|72|) is not independent and can be derived from the Einstein equations. 

In our model, as in [22j, we consider that the brane universe undergoes a natural inflation. The scalar field with 
an effective potential V(ip) appears in the bulk on account of the dynamics of the bulk gravitational field. As the 
universe inflates, the potential decreases and finally becomes negligible, so that the spatially flat universe is recovered 
in the low energy limit, with the bulk cosmological constant settling down to a stable value. However, during the 
inflationary phase, the bulk cosmological constant was presumably time-varying. All these considerations make the 
problem very hard to tackle, especially by analytical methods, since neither the potential, nor the bulk cosmological 
constant, can be tuned to zero. Attempts have been made to address these type of complicated problems with the 
help of numerical methods. In spite of that we choose to proceed analytically and sec how much information we can 
extract from this analysis. 

1. Solutions 

Assuming a solution of the type (|5"5)) . and using the relations between the coefficients of -j^ on both sides of the 
third Einstein equation ([751 leads us to the following possibility: 

^ - 2i\fif. (77) 
Using this result in the second Einstein equation (|74[) . we arrive at the relation 

f' = k 2 f> where k 2 = ~ + -^r, (78) 

so that 
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Therefore the form of ip(t, y) in general will be given by 



V>(*,2/) = ^-k 2 t - ^? I k^-dt + 2iV3 f ^-dy + constant. (80) 
2i 2i J A at J ay 

If A(t, y) is of the type ([60]), then 

?/) = ^-{k 2 t - lnAi) + 2«a/3 / ^-dy + constant. (81) 
2i 7 dy 

Analysing the above equations we find that the variation of the scalar field along the extra dimension depends on the 
variation of the warping function along the same direction. Further, the time-variation of the scalar field is related 
to the time- variation of the extra-dimensional scale factor. It appears that the time- variation of the warping function 
may be related to the time- variation of the extra-dimensional scale factor in most cases. This picture is also physically 
consistent, since these time- variations are expected to dominate only during the early phases of evolution, for which 
these must be functions decaying with time. 

We must remember that, along with this we have a time- varying bulk cosmological constant and an effective 
potential V(ip), which also decays with time. Since ip — ip\ + ip%, we can write V(ip) — V(ipi) + V(ip2)- At points 
where the space-variation of the potential is negligible, viz. at sufficient distances from the source, we can neglect the 
ip2 part and assume that V(ip) = V(ipi). However, close to the source, this assumption is invalid during the earlier 
epochs of time, when neither ip nor V(ip) is negligible. As a toy model we therefore consider a given form of potential 
V(ip), like the sine-Gordon potential, which dynamically generates the brane universe in the form of a thick brane in 
the background warped geometry from nothing. The universe inflates asymtotically at the Hubble rate H{1) 1 where, 
I is the AdS radius 

l = \J- T (82) 

The associated brane tension a will be of the type 

a = a c T(t,y), (83) 

where a c is the critical tension that reproduces the original RS braneworld. We assume that the bulk potential varies 
very slowly in space and time, so that the braneworld is created automatically from the dynamics of the space-time. 
The sine-Gordon potential is given by 

V(ip) = Bj +5iCos^. (84) 

Here B 2 must be nonzero. For simplicity, we choose #2 = 1. Thus B\ gives the value of potential when ip is zero i.e. 
V(0) = 2B\. Under this condition we have tpi — — ■02- We consider the situation where the scalar field ip is close to 
zero. In the lowest order, we can put f/ 1 = 0, so that the bulk is AdS with an effective cosmological constant A e j^, 
which is given by 

A e// = A + 2k 2 5i, (85) 

from which we have 

B, = (86) 
2k z 

with |A| > 2k 2 J5i, so that A. e ff < 0, thereby preserving the AdS nature of the bulk. The effective AdS radius is now 

(87) 

with h 2 = 1*1. 

To find the expression for the bulk cosmological constant, we subtract ([73| from (f76]) and arrive at the relation 
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where H = J|. Assuming 2/ = ^, we get 

2A = i (2 (A - 2Hf -f 2 + 2/) - . (89) 

Therefore the AdS radius and the bulk cosmological constant are both time-dependent and also dependent on the 
extra-dimensional coordinate. 



2. Cosmological interpretations: 



Although the bulk spacetime ([5]) is not conformally flat, its 4D part is found to possess a vanishing Weyl tensor 
and therefore is conformally flat. The Einstein tensor of the bulk is diagonal only under the condition (|3"5)) , but it 
is always diagonal for its corresponding 4D part. In general, the Ricciscalar for both the bulk and the 4D part are 
negative, indicating that both are AdS, except under some special circumstances. The constraint for the 4D part to 
be AdS is obtained as follows: 



f 2 + 2H 2 > \f + 3Hf + H\ 
while for the full 5D spacetime, we must have 



3(/ 2 + 2H 2 ) > 



3[f + 3Hf + H 



HA 
~A~ 



2fA A 
A A 



and 



5/' ; 



> 



r- 



f'A> 



A 



(90) 

(91) 
(92) 



The bulk pressure is isotropic along the brane, but is different along the extra dimension. However, the 4D part has 
isotropic pressure. The 4D part of this 5D bulk represents a universe with pure radiation if 



Rf + Rf + R = 0. 
In that case R~ V~t, and (JM]) reduces to 4 ft 2 + 2ft - 1 = 0. 



(93) 



VI. EFFECT OF PLACING THE BRANE AND LOCALIZATION OF FIELDS 

Having determined the bulk solutions, we now consider the effect of placing the single 3-brane in the bulk spacetimes 
considered above. This amounts to evaluating the brane tension a in terms of the parameters appearing in the 
expression for the bulk fields. For that we need to introduce the contribution of the 3-brane into the bulk action, as 
follows: 

S b = [ d 5 xVSV b W S -^^-. (94) 

The delta function appearing above will be peaked at the location of the brane (y = yo). The Einstein equations and 
hence the corresponding energy-momentum tensor for the brane-bulk system will get modified due to the appearance 
of the delta function term a5(y). Integrating the scalar field equation for the massless scalar at the location of the 
brane and using the Einstein equations, along with the constraint tp' 2 — — -§cjr, we obtain 

a ~ a c . (95) 

However, for an inflating brane, the departure of a from a c is not negli gible [33]. 

We now briefly discuss the nature of localization of the fields [U>ll2l [2J, [5(| for the spacetimes considered above. 
To explore the issue of localization of the various types of matter fields, the standard model fields are assumed to 
be distributed evenly over the entire bulk. The behavior of the KK modes along the extra dimension then indicate 
whether the fields are localized on the brane or not. Although the scalar (spin 0) and gravity (spin 2) fields can be 
localized on the positive tension brane in the RS background, it is generally not the case for the massless fermions. 
The massless chiral fermions can be localized on the positive tension RS brane in presence of Yukawa couplings with 
the bulk scalar. With an increasing warp factor, both the chiral zero modes may be confined to the brane, even in 
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the absence of any coupling with the bulk field. However, all the standard model fields may not be localizable on a 
single brane through gravitational interactions only, for other types of scalar fields and warp factors. 

To check for graviton localization, the spectrum of linearized tensor perturbations (restricted to 4D around the 
classical solutions of the warped spacetime), is analysed. The transverse traceless (TT) modes of the gravitational 
fluctuations represent the gravitons on the brane. These zero modes are analysed under the RS gauge conditions. 
In the typical RS background, the graviton zero mode is found to be localised and normalisable. In the braneworld 
models with massless bulk scalars and having a sine-Gordon potential, localized zero modes appear on the brane. 
With the flat 3-brane of the original RS setup being replaced by a dynamical brane in the generalized background as 
the one considered by us, the zero-mode graviton fluctuation is not guaranteed to be localized on the brane for an 
AdS bulk. The massless graviton wave function turns out to be non-normalizable once the Z 2 symmetry is relaxed, 
and the 4D Newton's constant may not be finite. The choice of bulk potential, the nature of coupling and the type 
of warp factor is crucial for achieving a normalizable zero mode. 

VII. SUMMARY AND DISCUSSIONS 

In this paper, we have considered five-dimensional warp product space-times with a time-dependent warp factor 
and a non-compact fifth dimension. The warp factor reflects the confining role of the bulk cosmological constant to 
localize gravity at the brane through the curvature of the bulk. Since this process of localization may include some 
time-dependence during a particular epoch of time, we have considered a generalized of warp factor in the form of 
e 2 /(*iJ/) ) which depends both on time as well as on the extra-dimensional coordinate. The extra-dimensional scale 
factor is also a function of time and of the extra coordinate, so that we are dealing with a very general type of bulk. 
The energy eigen values are evaluated and the condition for the WEC and NEC to hold in the bulk are discussed. 
Coupling these conditions with the conservation equations, we have been able to comment on the nature of matter- 
energy tensor for this bulk and the type of warping function. Assuming a definite correlation between the the rate of 
variation of the warping function and the normalized rate of variation of the extra-dimensional scale factor, we have 
found that the warping function turns out to be of the product type for a radiative brane. 

At low energies, when gravity and particles remain confined to the brane, we have found that the extra-dimensional 
scale factor depends on the variation of the warping function along the extra dimension. In general, the bulk is not 
conformally flat, although it can be so if the extra-dimensional scale factor is only a function of time. For a particular 
choice of the warping function such that the 5D metric may be conformally flat, the Einstein tensor for the bulk and 
its 4D part are diagonal and the respective Ricciscalars are generally negative. The braneworld represents an universe 
in an epoch of uniform acceleration. For a stabilized bulk of constant curvature, our model represents the usual RS 
scenario. Finally, at high energies, we have considered the most general warped metric for the bulk to be sourced 
by scalar fields living in the bulk. For an AdS bulk sourced by a non-self-interacting massless scalar field, we have 
determined the solutions for the scalar field, the warp factor and the extra-dimensional scale factor. Next, considering 
a dilaton-like scalar field with a sine-Gordon type of potential, we have determined the effective cosmological constant 
of the AdS bulk in terms of the bulk potential and have found the solution for the scalar field. The background 
geometry obtained in this case is described by a negative, non-constant Ricciscalar. The energy-momentum tensor, 
arising from this geometry may violate the strong energy condition depending on the nature of warping, as well as on 
the effect of the extra dimension, but the weak and the null energy conditions will be satisfied. The bulk in general 
will be anisotropic, but the corresponding four-dimensional part is isotropic. Although the brane tension is close to 
the critical value which reproduces the original RS braneworlds, yet the difference is not negligible. For this inflating 
braneworld model, the zero mode graviton fluctuation may not be localized on the brane. 

A detailed analysis of the localization of gravitons and other gauge fields will be presented in a separate report, 
along with the corresponding stability analysis, the nature of bulk singularities and other related aspects. 

Acknowledgments 

The major part of this work was done in IUCAA, India under the associateship programme. SG gratefully acknowl- 
edges the warm hospitality and the facilities of work at IUCAA. 



[1] K. Akama, Lect. Notes Phys. 176, 267 (1982) [arXiv:hep-th/0001U3] 
[2] V. A. Rubakov and M. R. Shaposhnikov, Phys. Lett. B 125, 136 (1983) 



15 



[3] M. Green, J. Schwarz and E. Witten, Superstring Theory, Vol. I and II (Cambridge University Press, Cambridge, 1987) 
[4] J. Polchinski, String Theory, Vol. I and II, (Cambridge University Press, Cambridge, 1998) 
[5] P. Horava and E. Witten, Nucl. Phys. B 460, 506 (1996); 475, 94 (1996) 

[6] A. Lukas, B. A. Ovrut, K. S. Stelle, and D. Waldram, Phys. Rev. D 59, 086001 (1999); Nucl. Phys. B 552, 246 (1999) 
[7] R. Maartens and K. Koyama, Living Rev. Relativity 13, 5 (2010); larXiv: 1004. 3962 Iiep-th] 

[8] Brane World: New Perspectives in Cosmology, Ed. K. Maeda and M. Sasaki [Prog. Theor. Phys. Suppl. 148, 181 (2002)] 
[9] P. D. Mannheim, Brane- Localized Gravity (Singapore: World Scientific, 2005) 
[10] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 3370 (1999) ; ibid. 4690 

[11] N. Arkani-Hamed, S. Dimopoulos and D. Dvali, Phys. Lett. B 429, 263 (1998); Phys. Rev. D 59, 0860 (1999) 
G. Dvali and M. Shifman, Phys. Lett. B 396, 64 (1997) 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B 436, 257 (1998) 
Z. Kakushadze and S.-H.H. Tye, Nucl. Phys. B 548, 180 (1999) 

B. Bajc and G. Gabadadze, Phys. Lett. B 474, 282 (2000) 
S.B. Giddings, E. Katz and L. Randall, JHEP 0003 (2000) 023 

S. L. Dubovsky, V. A. Rubakov and P. G. Tinyakov, Phys. Rev. D 62, 105011 (2000) 

A. Karch and L. Randall, JHEP 0105 (2001) 008 

E. Kh. Akhmedov, Phys. Lett. B 521, 79 (2001) 

A. Kehagias and K. Tamvakis, Phys. Lett. B 504, 38 (2001) 

1. Brevik, K. Ghoroku, S.D. Odintsov and M. Yahiro, Phys. Rev. D 66, 064016 (2002) 
M. Gogberashvili and P. Midodashvili, Europhys. Lett. 61, 308 (2003) 

G. Gibbons, K. Maeda and Y. Takamizu, Phys. Lett. B 647, 1 (2007) 
R. Davies and D. P. George, Phys. Rev. D 76, 104010 (2007) 
Y. X. Liu, L. Zhao and Y. S. Duan, JHEP 0704, 097 (2007) 

N. Barbosa-Cendejas, A. Herrera-Aguilar, M. A. Reyes Santos and C. Schubert, Phys. Rev. D 77, 126013 (2008) 
Y. Liu, X. Zhang, L. Zhang and Y. Duan, JHEP 0802(2008)067 

R. Guerrero, A. Melfo, N. Pantoja and R. O. Rodriguez, Phys. Rev. D 81, 086004 (2010) 
L. J. S. Sousa, W. T. Cruz and C. A. S. Almeida, Phys. Lett. B 711, 97 (2012) 
[12] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 (2000) 

[13] A. Lukas, B. A. Ovrut, and D. Waldram, Phys. Rev. D 60, 086001 (1999); 61, 023506 (1999) 
[14] T. Barreiro and B. de Carlos, J. High Energy Phys. 03, 020 (2000) 

M. Brandle, A. Lukas, and B. A. Ovrut, Phys. Rev. D 63, 026003 (2001) 

J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64 123522 (2001) 
[15] P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 565 269 (2000); P. Binetruy, C. Deffayet, U. Ellwanger and D. 

Langlois, Phys. Lett. B 477 285 (2000) 
[16] J. M. Cline, C. Grojean, and G. Servant, Phys. Rev. Lett. 83 4245 (1999); J. M. Cline and J. Vinet, J. High Energy Phys. 

02, 042 (2002) 

[17] P. J. Steinhardt, Phys. Lett. B 462, 41 (1999) 

C. Csaki, M. Graesser, L. Randall, and J. Terning, Phys. Rev. D 62, 045015 (2000) 

[18] E. E. Flanagan, S. H. H. Tye, and I. Wasserman, Phys. Rev. D 62, 044039 (2000); H. Stoica, S. H. H. Tye, and I. 

Wasserman, Phys. Lett. B 482, 205 (2000) 
[19] R. Maartens, D. Wands, B. A. Bassett, and I. Heard, Phys. Rev. D 62, 041301 (2000) 

E. J. Copeland, A. R. Liddle, and J. E. Lidsey, Phys. Rev. D 64, 023509 (2001) 

P. R. Ashcroft, C. van de Bruck, and A. C. Davis, Phys. Rev. D 66, 121302 (2002) 
[20] P. Brax, C. van de Bruck, A. C. Davis and C. S. Rhodes, Phys. Rev. D, 67, 023512 (2003) 

T. Damour and K. Nordtvedt, Phys. Rev. Lett. 70, 2217 (1993) 

D. I. Santiago, D. Kalligas, and R. Wagoner, Phys. Rev. D 58, 124005 (1998) 
[21] T. Barreiro, B. de Carlos, and E. J. Copeland, Phys. Rev. D 58, 083513 (1998) 

K. Choi, H.B. Kim, and H.D. Kim, Mod. Phys. Lett. A 14, 125 (1999) 
G. Huey, P.J. Steinhardt, B.A. Ovrut, and D. Waldram, Phys. Lett. B 476, 379 (2000) 
[22] Y. Himemoto and M. Sasaki, Phys. Rev. D 63, 044015 (2001) 

[23] W. D. Goldberger and M. B. Wise, Phys. Rev. Lett. 83, 4922 (1999); Phys. Rev. D 60 107505 (1999) 
[24] T. Tanaka and X. Montes, Nucl. Phys. B 582, 259 (2000) 

P. Kanti, S. Lee and K. A. Olive, Phys. Rev. D 67, 024037 (2003) 

R. Koley and S. Kar, Class. Quantum Grav. 22, 753 (2005); ibid. 24, 79 (2007); Phys. Lett. B 623, 244 (2005) 

S. Pal and S. Kar, Gen. Relativ. Gravit. 41, 1165 (2009) 
[25] T. Chiba, Phys. Rev. D 62, 021502 (2000) 

T. Wiseman, Class. Quantum Grav. 19, 3083 (2002) 

J. Khoury and R.-J. Zhang, Phys. Rev. Lett. 89, 061302 (2002) 

S. Kanno and J. Soda, Phys. Rev. D 66, 043526 (2002) 
[26] S. Mukohyama and L. Kofman, Phys. Rev. D 65, 124025 (2002) 

[27] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 484, 112 (2000); 485 208 (2000) 

C. Deffayet, Phys. Lett. B 502, 199 (2001) 
[28] R. Gregory, V. A. Rubakov, S. M. Sibiryakov, Phys. Rev. Lett. 84, 5928 (2000) 

C. Charmousis, R. Gregory and V. A. Rubakov, Phys. Rev. D 62, 067505 (2000) 



16 



P. Kraus, J. High Eneryg Phys. 12, Oil (1999) 
D. Ida, J. High Energy Phys. 09, 014 (2000) 

H. Kodama, A. Ishibashi, and O. Seto, Phys. Rev. D 62, 064022 (2000) 

C. V. D. Br uck, M. Dorca, R . Brandenberger, and A. Lukas, Phys. Rev. D 62, 123515 (2000) 
Y. Shtanov, |hep-th/0005193| 

J. Garriga and M. Sasaki, Phys. Rev. D 62, 043523 (2000) 
[29] L. Anchordoqui, C. Nunez, and K. Olsen, J. High Energy Phys. 10, 050 (2000) 
H. A. Chamblin and H. S. Reall, Nucl. Phys. B 562, 133 (1999) 
N. Kaloper, Phys. Rev. D 60, 123506 (1999) 
T. Nihei, Phys. Lett. B 465, 81 (1999) 

H. B. Kim and H. D. Kim, Phys. Rev. D 61, 064003 (2000) 

S. W. Hawking, T. Hertog, and H. S. Reall, Phys. Rev. D 62, 043501 (2000) 

S. Nojiri and D. Odintsov, Phys. Lett. B 484, 119 (2000) 
[30] F. Bonjour, C. Charmousis and R. Gregory, Class. Quantum Grav. 16, 2427 (1999) 
[31] J. Garriga, O. Pujolas and T. Tanaka, Nucl. Phys. B 605 192 (2001) 
[32] J. Garriga, O. Pujolas and T. Tanaka, Nucl. Phys. B 655 127 (2003) 
[33] K. Koyama and K. Takahashi, Phys. Rev. D 68, 103512 (2003) 

[34] P. Kanti, I. I. Kogan, K. A. Olive, and M. Pospelov, Phys. Rev. D 61, 106004 (2000) 
[35] T. Shiromizu, K-I. Maeda and M. Sasaki, Phys. Rev. D 62, 024012 (2000) 
R. Maartens, Phys. Rev. D 62, 084023 (2000) 

D. Ida, J. High Energy Phys. JHEP 0009 (2000) 014 

[36] P. Bowcock, C. Charmousis and R. Gregory, Class. Quantum Grav. 17, 4745 (2000); R. Gregory and J. A. Harvey, Class. 
Quantum Grav. 20, L231 (2003) 
M. Giovannini, Phys. Rev. D 76, 124017 (2007) 

G. Niz, A. Padilla and H. K. Kunduri, J. Cosmol. Astropart. Phys. 0804 (2008) 012 
[37] P. Kanti, S. Lee and K. A. Olive, Phys. Rev. D 67, 024037 (2003) 
[38] P. K. Townsend and M. N. R. Wohlfarth, Phys. Rev. Lett. 91, 061302 (2003) 

N. Ohta, Phys. Rev. Lett. 91, 061303 (2003) 
[39] S. Guha and S. Chakraborty, Int. J. Theor. Phys. 51 55 (2012) 

ibid. Gen. Relativ. Gravit. 42, 1739 (2010) 
[40] S. Ghosh and S. Kar, Phys. Rev. D 80, 064024 (2009) 

S. Ghosh, A. Dasgupta and S. Kar, Phys. Rev. D 83, 084001 (2011) 
[41] M. Gogberashvili, A. Herrera-Aguilar and D. Malagon-Morejon, Class. Quantum Grav. 29, 025007 (2012) 
[42] R. Kostera and M. Postma, JCAP 12, 015 (2011) 

I. P. Neupane, Phys. Rev. Lett. 98, 061301 (2007); Nucl. Phys. B 847, 549 (2011) 

[43] G. W. Gibbons, H. Lu and C. N. Pope, Phys. Rev. Lett. 94, 131602 (2005); W. Chen, Z. W. Chong, G. W. Gibbons, H. 
Lu and C. N. Pope, Nucl. Phys. B 732, 118 (2006) 

K. Maeda, N. Ohta, M. Tanabe and R. Wakebe, JHEP 0906 (2009) 036 
P. Binetruy, M. Sasaki and K. Uzawa, Phys. Rev. D 80, 026001 (2009) 
G. W. Gibbons and K. Maeda, Phys. Rev. Lett. 104, 131101 (2010) 

M. Minamitsuji and K. Uzawa, Phys. Rev. D 83, 086002 (2011); 84, 126006 (2011); K. Maeda and K. Uzawa, Phys. Rev. 

D 85, 086004 (2012) 
[44] A. K. Raychaudhuri and G. Mukherjee, Phys. Rev. Lett. 59, 1504 (1987) 

V. Dzhunushaliev, V. Folomeev and M. Minamitsuji, Rep. Prog. Phys. 73, 066901 (2010) 

D. P. George, M. Trodden and R. R. Volkas, J. High Energy Phys. 0902(2009)035 
[45] O. DeWolfe, D.Z. Freedman, S.S. Gubser and A. Karch, Phys. Rev. D 62, 046008 (2000) 
[46] A. Lue, R. Scoccimarro and G. D. Starkmann, Phys. Rev. D 69, 124015 (2004) 
[47] I. P. Neupane, Phys. Rev. D 83, 086004 (2011); Class. Quantum Grav. 26, 195008 (2009) 
[48] T. Gherghetta and A. Pomarol, Nucl.Phys. B 586, 141 (2000) 

R. Altendorfer, J. Bagger and D. Nemeschansky, Phys. Rev. D 63, 125025 (2001) 

A. Falkowski, Z. Lalak and S. Pokorski, Phys. Lett. B 491, 172 (2000) 
[49] K. Maeda and D. Wands, Phys. Rev. D 62, 124009 (2000) 
[50] C. Ringeval, P. Peter and J.-P. Uzan, Phys. Rev. D 65, 044016 (2002) 
[51] J. Carot and J. da Costa, Class. Quantum Grav. 10, 461 (1993) 
[52] P. Kanti, I. Kogan, K. A. Olive and M. Pospelov, Phys. Lett. B 468, 31 (1999) 

[53] M. D. Maia, E. M. Monte, J. M. F. Maia and J. S. Alcaniz, Class. Quantum Grav. 22, 1623 (2005); M. D. Maia, E. M. 

Monte and J. M. F. Maia, Phys. Lett. B 585, 11 (2004) 
[54] R. M. Wald, General Relativity (University of Chicago Press, Chicago, 1984) 

[55] G. S. Hall, J. Phys. A, 9, 541 (1976); C. A. Kolassis, N. O. Santos and D. Tsoubelis, Class. Quantum Grav. 5, 1329 (1988) 



